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SUMMARY ’ 


A review is given of recent advances in two distinct computational 
methods for evaluating turbulence fields, namely, statistical Reynolds stress 
modeling and turbulence simulation, where large eddies are followed in time. 
It is shown that evaluation of the mean Reynolds stresses, rather than use of 
a scalar eddy viscosity, permits an explanation of streamline curvature 
effects found in several experiments. Turbulence simulation, with a new vol- 
ume averaging technique and third-order accurate finite-difference computing 
is shown to predict the decay of isotropic turbulence in incompressible flow 
with rather modest computer storage requirements, even at Reynolds numbers 
of aerodynamic Interest. 


INTRODUCTION 


It is generally accepted that the Navler-Stokes equations represent the 
physics of Newtonian fluid flow fields at points in space and time away from 
sharp discontinuities such as shock waves. The intense shear layers present 
in turbulence fields, for example, the superlayer, extend over a much larger 
space than do shock waves, so that it can be expected that the Navler-Stokes 
equations will apply uniformly over the turbulence fields. Mathematically, 
turbulence can be considered, then, as the very complex solutions to these 
equations reflecting their nonlinear character and random sets of initial 
and boundary conditions. Because of the large set of possible solutions, 
unique solutions in turbulent flow can be expressed only in terms of averaged 
quantities over the time and/or space in the expectation that the turbulence 
is a slowly varying or stationary random process. 

The nonlinear character of the equations forces the use of numerical 
computations, yet attainment of these solutions of the time-dependent Navier- 
Stokes equations is not possible now or in the foreseeable future. The 
reason for this is contained in certain features of turbulence that are 
inimical to finite difference computations. At a point in space, turbulent 
flow is irregular in time, random in character. In addition, it is three- 
dimensional and composed of different structures having a large range of 
length and time scales. To compute the irregularities, it is necessary to 
employ computational methods that are accurately time dependent. Techniques 
appropriate to asymptotic solutions that permit overshoots or employ arti- 
ficial damping cannot be employed because they introduce errors in the 
Instantaneous flow-field development that are likely to grow. The 
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small-scale dimensions, also. Impose numerical stability or accuracy 
restraints on the allowable advancing time steps in the time-dependent solu- 
tions. For example," in an aerodynamic boundary layer, the smallest signifi- 
cant scale requires mesh spacing 10~^ of the boundary-layer thickness. The 
three-dimensionality of the flow, together with the small scale, requires on 
the order of 10l7 mesh points to define the flow field of an entire aircraft. 
The. corresponding allowable time step is 1 microsecond of real time. 

Clearly, these requirements for storage capacity and computer speeds are 
many orders of magnitude bhyond the best of current computers or those con- 
templated in the foreseeable future. To circumvent these difficulties it is 
necessary to resort to "turbulence modeling." 

The term "turbulence modeling" involves two distinct stages. The first 
eliminates the small-scale structure of the dependent variable through an 
averaging of some sort 3(ensemble, space, or time) that can be carried out 
with mathematical rigorj. The new averaged dependent variables are relatively 
slow varying so they can be resolved with finite difference methods utiliz- 
ing mesh dimensions and time intervals compatible with current computers. 

The nonlinear terms of the basic equations, however, introduce averages of 
moments of the dependent variables that must be expressed in terms of the 
averaged quantities retained to avoid a proliferation of dependent variables 
in excess of the numbers of equations available (the "closure" problem). The 
manner of doing this is the second stage of the "turbulence modeling" process 
and in this stage, considerable reliance is placed on comparisons with exper- 
imental data to compensate for the "physics" lost in the earlier averaging 
process . 

The name "Statistical Theory of Inhomogenous Turbulence," references 1 
and 2, has been applied to the process where averaging is performed at a 
point in space over a period of time long compared to the time scales of the 
largest eddies of the turbulence. The method is most suitable to steady- 
state mean flow fields, although it can be applied as well to flow fields 
that are slowly varying compared to the large eddy time scales. The result- 
ing dependent variables represent the slowly varying, mean flow field and are 
completely devoid of eddy structure. An alternative method, sometimes called 
"turbulence simulation," depends on averaging the Navier-Stokes equations 
over space volumes that are smaller than the largest eddies, but much larger 
than the tiniest eddies. The new dependent variables (the volume averages) 
retain their time dependence, but possess scales resolvable by current com- 
puter techniques. 

The NASA Ames Research Center has been sponsoring research in both of 
these areas, It is believed that in the short term most practical aerodyna- 
mic computations will be based on the statistical theory of turbulence. The 
aerodynamic 1st is interested in the mean properties of the flow fields with 
which he deals so that the behavior of individual eddies is more information 
than he needs. On the other hand, the modeling hypotheses made in the sta- 
tistical theory are largely Intuitive and supported rather tenuously by a 
meager supply of data. It hasn’t been possible to provide experimental 
Information on all of the multitude of variables that appear in the modeling 
equations. Turbulence simulation, once verified, will be able to be applied 
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to a variety of simple flow fields that emphasize particular modeling 
effects. For example, . simulation of homogeneous turbulence decay can shed . 
light on modeled turbulence dissipation in the statistical theory. Simple 
shear flows can yield information on the diffusion of turbulence kinetic 
energy, or of pressure fluctuations,. Thus, turbulence simulation will pro- 
vide an excellent basis for guiding statistical modeling. Ultimately, how- 
ever, the technique of turbulence simulation may go beyond merely providing 
modeling assistance to statistical theory and become the basis, for practical 
turbulence computation. Inherent in turbulence simulation is its direct 
evaluation of the large .eddies that are characteristic of particular flow 
fields and modeling of the small, scales that, are known to be universal in 
character. These are certainly the elements of a practical prediction 
scheme. 

In this .paper, a review will be given of some recent progress in the 
work sponsored by Ames. The recent contributions to the statistical theory ; 
of turbulence were performed at DCW Industries under Contract NAS2-8192. \ 

The work in turbulence simulation has been conducted at Stanford University 
under Grant NGR 05-020-622. 


. . SYMBOLS 

e,q2 specific kinetic energy of turbulence 

E power spectral density 

f function ■ , 

k wave number . 


L 

M 

P 



t,t* 


T 

u,v,w 

u, ■ 

1 

u 

T 

U 

o 


dimension of flow volume 
turbulence generation grid mesh spacing 
pressure ..... 

rate of strain of mean motion 

time, dummy variable of integration over time period 

period of integration > > turbulence scale 
velocity components in. x,y,z directions 

component of . velocity ini ith direction 

friction velocity, /r^/p 
.channel mean velocity 
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y 

r 

6 



A 

A 

avg 
e or 

V 

< / 
P 

T 


coordinate in ith direction 
distance normal to surface 
irrotational strain 
boundary- layer thickness , 

Kronecker 6 
mesh spacing 

dimension for space averaging (see table 1) , 

eddy diffusivity or eddy kinematic viscosity 
kinematic viscosity 
Von Karman constant 
density 

turbulent shear in boundary layer 

turbulent shear stress in direction 1 on surface normal to j 


T wall shear 

w 

X parameter in dissipation term of equation (16) 

0 ) Saffman's "pseudo vortlcity" or dissipation rate 

Superscripts; 

( ),< > average quantity 

( ) ' Instantaneous fluctuating quantity 


Subscripts: 

e edge of boundary layer 

l,j,k,m,n tensor indices 


STATISTICAL THEORY OF INHOMOGENOUS TURBULENCE 

To place the Wilcox contribution (ref. 3) into the proper context of 
current turbulence modeling it is necessary- to review briefly the basis of 
the statistical theory of turbulence.' For simplicity, attention will be con- 
fined to incompressible fluid flow. The Navier-Stokes equations are 
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where the notation of repeated indices represents summation over all the 
coordinate directions. The local velocity components, or pressure, are 
written as 

Ui(xj,t j 


which allows for slow variation of the mean quantity u^ in time t and 
expresses the fluctuating component within the integration period relative 
to the mean over the period. The" averaging process is defined as 


t+T/2 


dt 


t-T/2 

With equation (2), equation (3) becomes 
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which reduces to the requirement , , 
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When equation (1) is averaged according to equations (3) and (5), there 
results 
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Before equation (6) can be solved it is necessary to express the moments ' 

(the Reynolds stresses) in terms of the mean flow Uj and the characteristics 
of the turbul ence. If equation (1) is manipulated to yield a diffe rentia l 
equation for vij[uj , reference 1, in addition to, terms contai ning and 

Uj, new quantities involving quadra tures of uJ and p' and u|uju^ appear. 
An additional set of equations for Uj^uju^ leads to many new higher moments. 

In general, equations for higher moments lead to more new dependent variables 
than new equations. The means pf stopping this diverging process is achieved 
through "closure" hypothesis. 
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The most widely used closure method has been based on the concept of an 
eddy viscosity; for exam ple , where it is hypothesized that the turbulence 
Reynolds stresses u^Uj are proportional to the rate of strain of the mean 
flow, with the factor or proportionality, the eddy viscosity, being a scalar 
quantity. The similarity to Stokes stress relations in laminar flow and the 
role of viscosity is evident. With the eddy viscosity, e. 


ulul = 2 eS.. - 2/3 6. .e 
13 ij ij 


where the mean rates of strain, S.., are 
■ ij 


1 / ^'^l ■ 

^ij 2 ^ 8x^ 


(7) 


( 8 ) 


for ah incompressible ^fluid. Turbulence kinetic energy, e, is 


1 -r-r 
e = 2 u .u^ 


(9) 


The earliest theories, due to Prandtl, expressed 

e = E^X2,S^j ,i<j (iO) 

through an algebraic expression with k an, empirical constant.. Since this 
"closed" the problem at the level of equation (6) , it has been termed "first- 
order closure.". The method was quite successful for pipe flows and boundary 
layers where changes occurred gradually. In fact, extensions of this method 
to compressibility has provided several decades, of usefulness through accu- 
rate prediction of skin friction at high tfach numbers for two-dimensional 
aerodynamic configurations, evaluations of aerodynamic heating phenomena, 
and the assessment of concepts such as ■ transpiration cooling or ablation for 
alleviating the effects of aerodynamic heating. The method is most accurate 
when restricted to simple shapes. Its cost effectiveness in numerical com- 
putations, however, has made its use attractive in heuristic investigations 
even for conditions where.it cannot be expected to be accurate, namely, the 
zones of interaction between a shock wave and a turbulent boundary layer 
that is the topic of several of the papers in this conference. 


For. more general body shapes , especially those, having rapidly varying 
boundary- layer edge conditions, the introduction of rate processes in the 
development in the. turbulence is necessary. Here the eddy viscosity is 
expressed as a function of two turbulence quantities equivalent to velocity 
and time or length scales. As an example, the Saffman model (ref. 4) sets 

_ e^ _ turbulent specific kinetic energy . . 

^ 0 ) "pseudo vorticity" or "dissipation rate" ' ■ 
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where 


"and 


De 

Dt 


Ei(e,»,Sy) 


( 12 ) 



(13) 


The functions and f_ are not the exact terms that can be derived from 

equation (1), but are "modeled" to express the physical quantities actually 
present, rates of turbulence production, dissipation, and diffusion, in func- 
tional forms containing a closed set of dependent variables. These funcr 
tlonal forms are dimensionally consistent and contain .empirical constants 
that are expected to be insensitive to the character of Individual flow 
fields. A review of a large number of such two-equation models, where the 
modeling takes place in the equations for the turbulence quantities and is 
called "second-order closure," is given in reference 5. 


There are flow conditions, especially cases of three-dimensional flows, 
where the concept of the scalar eddy viscosity fai ls. F or these conditions 
it is necessary to evaluate the, Reynolds stresses directly. Several 

investigators have devised models of the direct Reynolds stress equations, 
references 6 through 8. The method described in this paper possesses , 

desirable features. It is relatively simple, reducing to the simple Saffman 
eddy viscosity model when flow conditions permit, and It Is easily extendable 
to compressible flows through the use of mass averaged dependent variables, 
reference 9. Furthermore, as Is shown later. It has proved to be quite suc- 
cessful In explaining the streamline curvature effects documented In Bradshaw's 
monograph, reference 10. 

.' ■ ■ ■ i 

The equation representing the Reynolds stresses of equation (6) can be 
written exactly as 
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(Molecular and Turbulent Diffusion of Stress) 
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I 

The terms on the left represent the convection of the Reynolds stresses. The 
first two terms on the right are the production terms; for example, where the 
mean flow strain interacts with the Reynolds stresses. The third term on the 
right is the dissipation term. It depends on quadratures of the instantaneous 
velocity gradients and therefore is largely dependent on the small scales of 
the turbulence. The next term, involving correlations of pressure and 
velocity gradient correlations, permits exchange between the individual 
Reynolds stresses and is termed the "tendency toward isotropy." This term 
vanishes for incompressible flow when equation (14) is contracted to yield 
an equation for the turbulent kinetic energy. The next term is the diffusion 
of pressure fluctuations by turbulence. Finally, the last term on the right, 
a divergence form, represents the diffusion of the Reynolds stresses by 
molecular and turbulent mechanisms. Dave Wilcox of DOT Industries models the 
Reynolds stress equation as 
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The turbulent dissipation rate, w, is given by 
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after Saffman, reference A. The turbulent kinetic energy, e, obtained as 
the contraction of equation (16) , is 

If- + - B*u)le + [t - 2 - S 1 

3t K 9x, L mnmn J [_ij wijj 3x. 


8*k 




(1) . 


(18) 


where 


a = 0.2638 
3 = 0.18 
a = 0.50 


3* = a*^ = 8/81 = .0988 


o = 0.5 


(9/2) 3* 
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These constants are consistent with the Saffman, formulation. Wilcox Intro- 
'duced an additional constant A*, which he evaluated by requiring the normal 

stresses in a flat plate boundary layer to have the ratio u'^ : v'^ : w'^ 
equal to 4:2:3 to conform to experiment. 

A comparison of the exact equation (14) and the modeled equation (16) 
indicates that modeling has been applied to all the terms other than turbu- 
lence production. The form of the modeled terms are expressed in a limited 
invariant tensor form, as compared to reference 7. Consistent with the 
Saffman model, the dissipation term does not represent the energy dissipated 
by molecular processes at the smallest scales, but rather accounts for the 
energy cascaded tq the smallest eddies by the larger eddies. Saffman.^ 
assumed that the scales of the smallest eddies adjust to accoimpdate the 
turbulence passed to them, and thereby avoids the need for two length 
scales. The quantity 1 - x + in equation (16) is always greater than 
zero, giving dissipation its proper sign. One can consider that the terms 
labeled "tendency to isotropy" and "production redistribution" together 
model the "tendency to isotropy" term of equation (14) as the contraction of 
all these terms is zero. The current Wilcox model neglects the "turbulent 
^diffusion of pressure." Finally, the "turbulent diffusion of stress," the 
; third order correlation, is modeled through an eddy viscosity assumption 
supplemented by an effective turbulent Prandtl number, a*. . 

One of the first problems to which this model was applied was the 
effect of streamline curvature in attached turbulent boundary layers. Cbn- 
aidefed was a constant pressure flow over a curved surface having a radius 
of curvature R in the plane of the two-dimensional mean flow. / From a 
perturbation solution of equations (6) and (16) to (18) corresponding to the 
wall region of the boundary layer, Wilcox found that mean velocity could be 
expressed in the usual wall layer dimensionless quantities as 


with 




Constant ' 


- (19) 


<[. = 15.4 
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The corresponding turbulence quantities were found to be 



t . ■ ■ ,'iV 
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In addition, the local shear turned but to be expressible in terms p£' an 
eddy viscosity ’ 


- 1 9u 

^ to 3y 


(26) 


These results are very significant In. several ways. First, the law of. the 
wall expression, equation (19) , has the same form found to correlate experi- 
mental data measured over a curved, wall. The value of found experiment- 
ally by R. N. Meroney, was. <|) = 12, which is in reasonable agreement ' with the 
modeled value of equation (20) . Furthermore, the specific kinetic energy of 
the turbulehcej equation (21), is shown to be relatively unaffected 'by the 
streamline. curvature. Again, this is in agreement with experimental data 
for corresponding flow conditions where So and Mellor, reference 11, found 
C mav / u^ unaffected by curvature for convex walls. A comparison of 

equations (23) through (25) shows that streamline curvature affects the 
individual normal stresses differently. The normal stresses along and 
across the surface increase, while the normal stress perpendicular to the 
surface decreases and in an amount larger than the others. The difference 
in sign of the behavior of the normal stresses demonstrates why the specific 
kinetic energy remains essentially constant. It is interesting that the 
expression equation (26) for the shear stress, t, shows an eddy viscosity 
form, but where the characteristic turbulence velocity scale is n ot t he 
specific kinetic_ energy, e, reference 4, but the normal stress v*^ . The 
sensitivity of v^ to the effects of streamline curvature and its direct 
effect on. the shear explains much of the streamline curvature effect dis- 
cussed by Bradshaw, reference 10. * 


Another comparison of the Wilcox Reynolds stress computations with 
experimental data is shown in figure 1. These data were obtained by 
•Bradshaw, reference 12, in. a turbulent boundary layer on a flat plate on 
which an adverse pressure was suddenly applied. The data Include profiles 
of the three individual normal Reynolds stresses, the shear stress, the 
specific kinetic energy, and the. mean velocity at a station 1.52 m after the 
adverse pressure gradient was applied. To initiate the computations, 
extensive use was made of the Bradshaw mean flow and hot-wire data at the 
start of the adverse pressure gradient. The initial mean velocity and 
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specific kinetic energy profiles were directly obtained from the data. The 
w profile was deduced from the measured' turbulent shear stress and mean 
velocity profiles by utilizing the eddy viscosity relationship 


e du 
^ ^ 0 ) 3y 


(27) 


Because the shear t and 3u/3y both approach zero near the boundary-layer 
edge, the determination of «j there becomes uncertain and is a source of 
error. The value of -oi found at the boundary layer edge and used at all 
stations as a boundary condition on equation (17) was Ug = 10“** Ug/v. 

The Reynolds stress theory, shown in figure 1 as the solid lines, -does, 
a good job in predicting the ma ximum lev els of the normal stresses, however, 
the shap es of the profiles for and w'^ are missed badly. In. contrast, 

the v'^, the normal stress perpendicular to the surface, is predicted quite 
well everywhere in the boundary layer.’ ^It should be no ted that af ter e xper- 
leticing the adverse pressure gradient*’ the ratio u ' 2^^ ; v' ^^ v : ^ 

9:5:7, which the theory seems to predict without modification of the con- 
stant X*. The comparisons with the data for the shear stress, turbulent 
energy, and mean velocity also contain the predictions by the eddy viscosity 
model of Saffman indicated by the dashed lines. Because the Wilcox model 
relies so heavily on the Saffman model, it is not surprising that the two 
theories . give such similar results*, especially at a station where the turbu- 
lence is nearly in equilibrium. For these quantities, the predictions, 
while encouraging, _ yield results that can stand improvement.- 

From this brief comparison with data, it can be concluded that the 
Wilcox Reynolds stress model with second-order closure predicts the effects 
of streamline curvature rather accurately in the wall region of a boundary 
layer. Of the normal stresses, the is most strongly affected. This, in 

turn, affects the shear stress since in this region T ~ 1/X*(v' ^/w) (3u/9y) . 
The latter is a most important practical, result in that it shows that- current 
eddy viscosity models with v'^ replacing the turbulence energy as the charac- 
teristic turbulence velocity still may apply in regions of streamline curva- 
ture.' Despite the success of the model near the surface, modeling improvements 
are required to increase the accuracy of predictions away from the surface. 
Perhaps the diffusion of pressure fluctuations by turbulence cannot be neg- 
lected. Further, it may be necessary to modify the constants of the Saffman 
theory; for particular flow fields to enhance accuracy at the sacrifice of 
generality.' ' . • . i . 

TURBULENCE SIMULATION 


The turbulence simulation calculations to be described here are being 
conducted at, Stanford University, with ■Profs. William Reynolds and Joel 
Ferziger as the principal Investigators. The primary goal of the Stanford 
activity is to derive turbulence simulation techniques ultimately suitable 
for use with compressible fluids and complex aerodynamic configurations. The 
initial stages of the program, however, have been confined to very simple 
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flow fields and to incompressible fluids so that comparisons of the tech- 
niques employed could be made with the results. of other workers in turbu- 
lence simulation and with certain fundamental experiments based on hot-wire 
measurements. 

The program began with a careful reconsideration of the volume averag- 
ing process by Leonard, reference 13. The new averaged equations contained 
second- order terms that required finite difference methods with truncation 
errors ..fourth order in space and second order in time. These equations have 
been used to study the decay of incompressible, isotropic flow with differ- 
ent numerical techniques and alternative models for subgrid closure. The 
work is being extended, currently, to flows with irrotational plane strain, 
simple shear, and a jet into- an axial stream. 

■ Table 1 and figure 2 compare the , bases of the averaging technique 
devised by Leonard with those of the .more • conventional approach. The 
averaged Navier-Stokes equation shown at the top of table 1 is common to ^ 
both methods. The dependent variables, however, have different' meaning. 

The conventional average quantity is ,.a uniformly weighted average over a... 
volume of dimension Aj around the point x j . The volume dimension A j • is 
customarily set equal to the finite difference mesh dimension.. The Leonard 
average is weighted toward the center of the averaging volume Xj as 
G(xj-Xj') is required to decay to a zero value with increasing distance 
xj-xj.'. Specifically, in most of the numerical work a Gaussian weighting 
function has been adopted with A^vg» not. necessarily set equal to the mesh, 
spacing, as a parameter. In either averaging method, the gradient of an 
average equals the average of a gradient. Key to the difference in the 

averaging techniques is the manner in which the local velocity vector is 
divided into the resolvable and irresolvable eddy contributions. This is 
shown mathematically in table 1 and diagramatically in the figure 2. In the 
conventional case, the resolvable velocity is identified with the mean velo- 
city In the box so that the irresolvable component , then, averages to zero. 
With this definition the averaged nonlinear moment can be expressed as the 
product of the individual average velocity components plus the Reynolds 
stress. It is emphasized that this Reynolds stress is time dependent, just 
as are the averaged velocities related to the large eddies. In Leonard 
averaging, the resolvable component of velocity u^ is treated as a con- 
tinuously varying quantity across the averaging volume, and the irresolvable 
fluctuations are defined relative to the local Uj^. The averaged nonlinear 
moment is much more complex containing a term representing the divergence of 
the quadrature of the mean velocities as well as a series of quantities that 
are identified in the aggregate as a Reynolds stress. The modeling of the 
Reynolds stresses, although they are inherently different quantities, is the 
same functionally in either method. The Reynolds stress of the unresolved 
components, Rij , for either averaging method is expressed in terms of an 
eddy viscosity and the instantaneous strain of the resolved eddies . The 
definition of the eddy viscosity follows Smagorinsky, reference 14, where the 
length scale is set equal to the mesh spacing and the frequency scale is 
related to the root of the sum of the squares of all the strain elements. 

An alternative eddy viscosity has been utilized by the Stanford group where- 
in the frequency scale depends on the .vorticity of the resolved eddies. 
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having the advantage of going to zero where the flow Is irrotatlonal. The 
latter approach may aid In defining the boundaries between turbulent and 
nonturbulent , irrotatlonal flow regions. . • . 

The main contribution resulting from the Leonard averaging process is 
the addition of the term A|^g/24(8^/9x?) (u^u.) to the expression for the 
averaged moment. Since this term is contained within an operator 9/9x| it 
is of third order when expressed in finite difference form. To retain the 
significance of this term, the Stanford group was forced to adopt finite dif- 
ference methods with fourth-order truncation errors in space,, reference .15 . . 
The usual second-order finite differencing in time was retained. 


The results of some. of the Stanford computations of the decay of homo- 
geneous, isotropic turbulence are shown in figures 3 through 6. The data used 
for comparison were measured by Compte-Bellot and Corrsin, ref erence -16 , in a 
uniform flow behind a turbulence generating grid at a Reynolds number per 
foot of- about 200,000. The experimental three-dimensional turbulent energy 
spectrum at a fixed position in the channel downstream of the grid, shown in 
figure 3, was adopted as the initial conditions for, the calculations. 

Because the mean flow rate is uniform,, the distance along the channel serves 
as the independent variable of the problem and can be expressed in time as . • 
tUp/M, where M is the mesh spacing of the turbulence generation grid . At 
the higher wave numbers, the spectrum possesses the k~5/3 distribution 
expected of homogeneous, isotropic turbulence according to the theory of 
Kolmogoroff, reference 17. The vertical lines in figure 3 define the range 
of wave numers, k, resolved in the calculations by mesh spacings where each, 
side of the calculational volume is divided into 16 or 32 parts ^ The lower. 


bound results from the imposition of periodic boundary conditions on opposite 
faces of the cubic volume,- so that 



k = — ' 
min L 


./ (28) 


where L is the length of the side of the cube. For the cube divided in 16^ 
parts, L ' is taken to be 24 cm. For the cube divided in 32^ parts, L is 
32 cm. The corresponding kjnin> fi’Pm equation (28), are 0.26 cm“^ and 0.2 
cm“l, respectively . The upper value of the resolved k is dependent on the 
computational mesh dimension, being smaller for the 32^ mesh division. 

Within these bounds, the 323 mesh division can resolve eddies that account 
for 80% of .-the total turbulence energy, whereas the 16 ^ mesh division only 
accounts for about 71% of the energy. - 


The next three figures (4 to 6) show the evolution of the turbulence 
spectra. at times, or distances beyond the turbulence generating grid, after 
the initial value set at the farthest upstream test station,. All of the fig- 
ures represent the case of the cube divided into 16^ parts and with computa- 
tional mesh spacings of 1.5 cm. The results found from the 323 case showed no 
marked imiprovements in the calculations.. The computation was initiated by 
expressing., the turbulence in a 3-dimenslonal Fourier series. The Fourier 
transform of this, velocity series to wave-number space then was matched 
directly to the experimental energy spectrum to give the magnitude of the 
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coefficients corresponding to each wave-number vector. These vectors were . 
then randomly oriented with a random number generator. Finally, an inverse 
Fourier transform was applied to iQchieve. the initial turbulence field compat- 
ible with the measured energy in physical space. The decay of the turbulence 
field was then computed using the time-dependent methods described .earlier.. 

Figure 4 shows the effect of setting Aayg = 0 utilizing ..the 

Smagorinsky subgrid eddy viscosity model (see table 1) in the. computations. 
The dashed. lines represent the data at the initial time, tUo/M = 42, and a 
subsequent time, tUp/M = 98. The computed spectra are represented with the 
solid lines. Initially the computations are made to agree with the 
tUo/M = 42 data line as described earlier. The development of the spectra 
after this time is the test of the method and the subgrid model. From table 
1, it is noted that with A^vg = 0, the weighting function in • the Leonard . . 
averaging is a Dirac function and the averaged Navler-Stokes equation is 
formally the same as in conventional-averaging, although the meaning of the 
dependent variable u^ is different. For Ag.^g = 0, as time progresses, a 
piling up of energy occurs at the high wave-number end that becomes prpgresr 
sively worse. - 

Figure 5 shows the evolution of the energy spectrum if the averaging 
volume dimension is Increased to equal to the computational mesh dimensions, 
Aavg ® A. Again, the Smagorinsky eddy viscosity model is used for the sub- 
grid Reynolds stresses. Note that the dashed lines representing the data • - 
are different from those in the previous figure because the data have been . 
averaged according to the Leonard fomula with A^^g ” ^ order to 

express the data in the same dependent variables as in the computations . 

The unusually high energy at the high wave numbers evident in the previous 
figure for “ 0 ^t tUp/M = 98 is reduced here for A^^g = 

Finally, in figure 6, the results of using A^^^^g = 2A- are shown. Again the 
data have been averaged to make the data and computation dependent variables 
correspond. The dashed line at tUg/M = 98 is the averaged data. The 
points in symbols show the computed spectrum at tUp/M = 98 for both the 
Smagorinsky and vorticity models of the subscale eddy viscosity (see' table ‘ 
1). There is little to choose between the different models for this 
example. Both show excellent agreement with the data, and the piling-up of 
the calculations at high wave number is no longer evident.' It appears that 
the use of most appropriate for this type of problem and that 

the additional term in the Leonard averaging process acts to redistribute 
the energy within the spectrum in an appropriate manner. 

j 

The last figure, figure 7, shows some early results of computations of 
the distortion of homogeneous turbulence by irrotational plane strain. A 
field of homogeneous isotropic, turbulence is pa:ssed through the transition 
section shown schematically in the upper left of the figure. , The cross- 
I sectional area is constant so that' w, the mean velocity ih‘the z direction 
remains constant. In the center portion, the walls perpendicular to the y- 
axls move in to just balance the outward movement of the walls perpendicular 
to the x-axis, resulting in a constant irrotational plane strain of the mean 
flow. .Downstream of the transition region, the walls again become parallel. 
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The development of the mean normal stresses with distance along the channel 
as computed by turbulence simulation is shown in the lower portion of the 
figure. Again, a cubic volume of 163 mesh points was employed. At station 
0 the turbulence, init i ated a s de scribed earlier for the decay problem, is 
nearly Isentroplc u' ^ = v'^ = w'^. In the region of con stant irrotatlorjal 
strain, t he n ormal stress in the expanding direction, u’^, diminishes, 
whereas, v’^, in the contracting direction increases. The axial normal 
stress remains relatively unchanged. In the downstream parallel section, the 
turbulence appears to be returning to Isotropy as expected, but very slowly. 

The results of turbulence simulation calculations described here are 
most encouraging. Through adjustment of the volume-averaging process. and 
the use of higher order finite differencing, it has been possible to accu- 
rately predict the decay rate and spectral content of the large eddies of 
homogeneous, isotropic turbulence in Incompressible flow. No difficulties, 
have, been encountered for the decaying type' of flow with the standard sub- 
grid modeling techniques. In addition, early attempts at handling more com- 
plex flows are yielding reasonable results with relatively little extra 
effort. Perhaps the most Important result to date is the attainment of 
accurate turbulence prediction at Reynolds numbers of aerodynamic Interest 
with the rather small number of mesh points employed. 


CONCLUDING REMARKS 


The success of the Reynolds stress modeling in the statistical theory 
of turbulence within a boundary layer in the immediate vicinity of the wall 
is most encouraging. Although improvements in modeling are still required 
for the wake region of the boundary layer, a way of doing this other than by 
reasoned intuition seems to be close at hand through turbulence simulation. 
In the near future, significant advances can be expected through the inter- 
section of these seemingly disparate approaches to the evaluation of turbu- 
lence flows. Solutions to simple problems at aerodynamic Reynolds numbers 
appear to be possible currently as computer storage requirements and subgrid 
modeling do not seem to be as critical as was anticipated earlier . 
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TABLE 1. COMPARISON OF CONVENTIONAL AND LEONARD AVERAGING 
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Figure 1. Comparison of Reynolds stress computations id.th experimental 

data of Bradshaw (adverse pressure) . 

















